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S '. ABSTRACT 

; I ' This paper gives a generalization of the AJL algorithm for quantum computation of the Jones polynomial to 
continuous ranges of values on the unit circle for the Jones parameter. We show that the Kauffman-Lomonaco 
3-strand algorithm for the Jones polynomial is a special case of this generalization of the AJL algorithm. 

^\ ' Keywords: knots, links, braids, quantum computing, unitary transformation, Jones polynomial, Temperley- 
CN , Lieb algebra 

H ■ 1- INTRODUCTION 

In^^ and in^** we gave a quantum algorithm for computing the Jones polynomial via a unitary representation 
, of the three-strand Artin braid group to the Temperley-Lieb algebra. In the bracket polynomial version of this 
representation (see Section 2 of the present paper) the representations were unitary for certain continuous ranges 
of choice of the polynomial variable A on the unit circle in the complex plane. In this paper we show that 
these three-strand representations are a subset of unitary representations of the Artin braid group on arbitrary 
numbers of strands and corresponding continuous ranges of the variable A on the unit circle. These more general 
representations are in fact generalizations of the AJL representations^'^ that were originally defined at certain 
, roots of unity in the unit circle. 

<N 

. The paper is organized as follows. In Section 2 we review the bracket polynomial model for the Jones 

' polynomial, and its relationship with representations of the Temperley-Leib algebra. In Section 3 we review the 3- 
. strand representation. In Section 4 we detail diagrammatically the construction of the generalized representation 
' and show how it is related to the 3-strand representation and to the AJL representation. In Section 5 we give 
. a diagrammatic proof of the requisite trace formula that is needed to make this representation into a quantum 
algorithm for computing the Jones polynomial. Much remains to be explored in these directions. The present 
paper was sparked by our work in^^ on NMR quantum computing, and there will be a sequel to the present 
paper^'' that relates the present work to NMR research. 



2. BRACKET AND TEMPERLEY LIEB ALGEBRA 

The bracket polynomial^ model for the Jones polynomial**^^' is usually described by the expansion 

(X) = A()<) + yi-i(X) (1) 

and we have 

{K O) = i-A' - A-'){K) (2) 

(X) = i-A'){-^) (3) 

{Y) = (4) 




Figure 1. Reidemeister Moves 



The bracket expansion of a knot or link diagram is invariant under Reidemeister moves II and III as shown 
in Figure 1, and can be normalized to be invariant under the first Reidemeister moves by multiplication by an 
appropriate power of —A^. Once normalized, it is a version of the Jones polynomial,^ differing from it by a 
simple change of variable. 

The key idea behind the present quantum algorithms to compute the Jones polynomial is to use unitary 
representations of the braid group derived from Temperley-Lieb algebra representations that take the form 



where CTj is a standard generator of the Artin braid group, ^ is a complex number of unit length, and Ui is a 
symmetric real matrix that is part of a representation of the Temperley-Lieb algebra. A diagrammatic version 
of the Temperley-Lieb algebra puts the form of this representation in exact correspondence with the bracket 



form correspond to the generator Ui of the algebra when y\ corresponds to the braid generator o-j. For 

more details about this strategy and the background information about the Jones polynomial, the bracket model 
for the Jones polynomial and the Tcmpcrlcy-Licb algebra the reader may wish to consult.^' In the 
following sections, we have made use of such diagrammatic techniques and have included some material to make 
the paper partly self-contained. 



3. TWO PROJECTORS AND A UNITARY REPRESENTATION OF THE THREE 

STRAND BRAID GROUP 

It is useful to think of the Temper ley Lieb algebra as generated by projections ej = Ui/d so that e| = ej and 
ejejiiej = rej where r = (5~^ and Cj and ej commute for \i — j\ > 1. 

With this in mind, consider elementary projectors e = and / = \B){B\. We assume that {A\A) = 

{B\B) = 1 so that 6^ = 6 and = /. Now note that 



p{a^) = AI + A-^Ui 




correspond to the identity element of the algebra and the arcs in the 




efe=\A){A\B){B\A){A\ 



{A\B){B\A)e = Te 



Thus 



e/e = re 



where t = {A\B){B\A). 



Further author information: L.H.K. E-mail: kaufrman@uic.edu, S.J.L. Jr.: E-mail: lomonaco@umbc.edu 



This algebra of two projectors is the simplest instance of a representation of the Temperley Lieb algebra. 
In particular, this means that a representation of the three-strand braid group is naturally associated with the 
algebra of two projectors. 



Quite specifically if we let (^| = (a, b) and \ A) = {a, by the transpose of this row vector, then 

e = \A){A\ = 

is a standard projector matrix when o,^ + 6^ = 1. To obtain a specific representation, let 
lei = 



ab 
ah 6^ 



" 1 
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and 62 = 


a? ah 








ab &\ 





ab 


and 6261 = 


■ a^ 


" 


6162 = 





ah 






It is easy to check that 616261 = a^6i and that 626162 = 0^62. 
Note also that 6162 

L u u 

We define Ui = Sci for i = 1, 2 with = S~'^. Then we have , for i = 1, 2 

Uf = 5Ui , U1U2U1 = f/i, U2U1U2 = U2. 

Thus we have a representation of the Temperley-Lieb algebra on three strands. See® for a discussion of the 
properties of the Temperley-Lieb algebra. 



Note also that we have 
while 



trace{Ui) = trace{U2) = S, 

trace{UiU2) = trace{U2Ui) = 



1 



where trace denotes the usual matrix trace. These formulas show that the trace is working correctly with respect 
to the bracket evaluation. See the last part of the present paper for a more extended discussion of this point. 

Now we return to the matrix parameters: Since + 6^ = 1 this means that + b^ = 1 whence 6^ = 1 — 6~^. 
Therefore h is real when 6"^ is greater than or equal to 1. 



We arc interested in the case where S 
stances the braid group representation 



-J^ — A ^ and A is a unit complex number. Under these circum- 
p{ai) = AI + A-^Ui 



will be unitary whenever Ui is a real symmetric matrix. Thus we will obtain a unitary representation of the 
three-strand braid group B3 when S"^ > 1. 

For any A with d = —A^ — A^"^ these formulas define a representation of the braid group. With A ~ exp{i9), 
we have d = —2cos{29). We find a specific range of angles in the following disjoint union of angular intervals 

e e [0, 7r/6] U [tt/S, 27r/3] U [57r/6, 77r/6] U [Aw/S, 57r/3] U [ll7r/6, 27r] 

that give unitary representations of the three-strand braid group. Thus a specialization of a more general repre- 
sention of the braid group gives rise to a continuous family of unitary representations of the braid group. 



3.1. A Quantum Algorithm for the Jones Polynomial on Three Strand Braids 



We gave above an example of a unitary representation of the three-strand braid group. In fact, we can use 
this representation to compute the Jones polynomial for closures of 3-braids, and therefore this representation 
provides a test case for the corresponding quantum computation. We now analyse this case by first making 
explicit how the bracket polynomial is computed from this representation. This unitary representation and its 
application to a quantum algorithm first appeard in.^^ When coupled with the Hadamard test, this algorithm 
gets values for the Jones polynomial in polynomial time in the same way as the AJL algorithm.^ It remains 
to be seen how fast these algorithms are in principle when asked to compute the polynomial itself rather than 
certain specializations of it. 



First recall that the representation depends on two matrices Ui and U2 with 



Ui 



S 




and U2 



6-6-^ 



The representation is given on the two braid generators by 



p{cri) = AI + A-^Ui 
p{(j2) = AI + A~^U2 



for any A with 5 = —A^ — A^^. and with A = exp{iO), then S = —2cos(26). We get the specific range of 
angles 6 £ [0, 7r/6] U [vr/S, 27r/3] U [5n/6, 7it/6]U [47r/3, 57r/3] U [ll7r/6, 2tt] that give unitary representations of the 
three-strand braid group. 



Note that tr{U\) = tr{U2) = S while tr{UiU2) = tr{U2Ui) = 1. If 6 is any braid, let I{b) denote the sum of 
the exponents in the braid word that expresses 6. For b a three-strand braid, it follows that 

p{b) = A'^^h + Tib) 

where / is the 2x2 identity matrix and rib) is a sum of products in the Temperley Lieb algebra involving Ui 
and U2- Since the Temperley Lieb algebra in this dimension is generated by U2, U1U2 and U2U1, it follows 
that 

(6) ^ A'^^^^5^ + tr{T{h)) 

where b denotes the standard braid closure of b, and the sharp brackets denote the bracket polynomial as described 
in previous sections. From this we see at once that 

(6) =tr(p(6))+A^W(<52_2). 

It follows from this calculation that the question of computing the bracket polynomial for the closure of the 
three-strand braid b is mathematically equivalent to the problem of computing the trace of the matrix p(b). 



The matrix in question is a product of unitary matrices, the quantum gates that we have associated with 
the braids cri and 02- The entries of the matrix p{b) are the results of preparation and detection for the two 
dimensional basis of qubits for our machine: 

mm- 

Given that the computer is prepared in the probability of observing it in state \i) is equal to | p. 
Thus we can, by running the quantum computation repeatedly, estimate the absolute squares of the entries of 
the matrix p{b). This will not yield the complex phase information that is needed for either the trace of the 
matrix or the absolute value of that trace. 



However, we do know how to write a quantura algorithm to compute the trace of a unitary matrix (via the 
Hadamard test). Since p{b) is unitary, we can use this approach to approximate the trace of p{b). This yields 
a quantum algorithim for the Jones polynomial for three-stand braids (evaluated at points A such that the 
representation is unitary). Knowing tr{p{b)) from the quantum computation, we then have the formula for the 
bracket, as above, 

(6) = trace{pib)) + A'^^^^S"^ - 2). 
Then the normalized polynomial, invariant under all three Reidemeister moves is given by 

/(6) = (-A3)-^W(6). 
Finally the Jones polynomial in its usual form is given by the formula 

vmt) = fmt-'/^). 

Thus we conclude that our quantum computer can approximate values of the Jones polynomial. 

4. GENERALIZING THE AJL REPRESENTATION 

In this section wc show how the KL (KaufFman-Lomonaco) algorithm described in the previous section becomes 
a special case of a generalization of the AJL algorithm: Here we use notation from the AJL paper. In that 
paper, the generators Ui (in our previous notation) for the Temperley-Lieb algebra, are denoted by Ei. We will 
first describe the AJL representation of the Tcmpcrley-Lieb Algebra and we will show how that representation 
works for a continuous range of values of the parameter 9 described below. In the original treatment of AJK^ the 
values are restricted to a discrete range corresponding to exp{i6/2) a root of unity. We observed this phenomena 
of continuous ranges of values in our work on the three-strand model, described in the previous section, and 
found a way to extend it to AJL, as will be detailed below. This work has benefited from interaction with all of 
the authors of the paper^^ and our collaboration in bringing these algorithms to application in NMR quantum 
computing. The contents of this section will be connected with NMR quantum computing in a paper that is 
under preparation.^^ 

In this section wc will construct matrix representations of the Temperley-Lieb algebra. In the next section, 
we will discuss the structure of trace functions on these representations that can be used to produce quantum 
algorithms. 

Let A/c = sin{k6). For the time being 9 is an arbitrary angle. Let A = iexp{i9/2) so that d = —A^ — A~^ = 
2cos{e). 

We need to choose 9 so that sin{k9) is non-negative for the range of fc's we use (these depend on the choice of 
line graph as in AJL). And we insist that sin{k9) is non-zero except for A; = 0. Then it follows from trigonometry 
that for all k 

d = (A/c_i -I- A/c+i)/ Afc. 

We shall see, below, that the values of k will range from to r for a fixed r > 2 in a given representation. We 
ask that sin{r9) be greater than zero and can take 9 in the continuous range < 6 < w/r. 

The AJL representation of the Tempcrlcy-Licb algebra is based on the complex vector space Hn whose basis 
is {|z)} where z is a (0, 1) bitstring of length n. Each bitstring i is seen as corresponding to a walk on a line-graph 
Gr with r — 2 edges and r — 1 nodes. For example, the graph below is G5. 

1 2 3 4 

Bitstrings represent walks on a line graph with corresponding to a step to the left and 1 corresponding to a 
step to the right. The walk begins at the left-most node, labelled 1. Thus 1011 represents the walk Right, Left, 
Right, Right ending at node number 3 in the line graph above. We say that a string of n bits is a walk on Gr if 
the walk remains inside Gr for all its steps (i.e. one never faces the instruction to go beyond the right-most or 
left-most end points of the graph). We let Hn^r denote the subspace of Hn spanned by bitstrings that are walks 
on Gr- 



The representation of the Temperley-Lieb algebra is denoted by 



$ : TLn — > Matr{Hn,r) 

where T_L„ denotes the n-strand Temperley-Lie algebra generated by /, • • • , En-i and Matr{Hn.r) denotes 
matrix mappings of this space in the given bitstring basis. The AJL representation is given in terms of Ei such 
that = dEi and the E^ satisfy the Temperley-Lieb relations. Each Ei acts non-trivially at the i and i + 1 
places in the bit-string basis for the space, lip denotes the bitstring and we are computing then we define 

z{i) to be the endpoint of the walk described by the bitstring p using only the first {i — 1) bits of p. Each Ei{p) 
is based upon \z{i)-i^^z{i)-:^z{i)+\- our example we have p = 1011 represents the walk Right, Left, Right, 
Right ending at node number 3 in 

1 2 3 4, 

and z{l) = 1, z{2) = 2, ziS) = 1, ^(4) = 1, ^(5) = 3. 

More precisely, if we let 

\v{a)) = i^Xa-l/Xa, ^\a+l/Xaf 

(i.e. this is a column vector. T denotes transpose.) Then 

E,^\v{zm{v{zm. 

Here it is understood that this refers to the action on the bitstrings 



and 

10 



obtained from the given bitstring by modifying the i and i + 1 places. The basis order is 01 before 10. 



The explicit form of the the transformations Ei is given by the equations below for a generic matrix E. We 
need to explicate ^{Ei) and we shall simply write Ei instead of ^{Ei). 



E = vv^ = 



A_/Ao ^ 

Ao 



Here = (a/t-) v y*^)- It is easy to see that 



A+/Ao 



E' 



> A_ + A_| 
^ Ao 



■)E 



since E'^ 



vv vv 



= {v'^v)E. For Ei, we take A- = Xz(i)-ii A+ — Xz(i)+i and Ao — \z{i)- For the action of 



Eiiji) for a given bitstring p, let p\i denote the restriction oip to the first i—1 bits in the string. Then we need to 
explicate Ei at the bitstring p and hence its values at p\i k I) where k and I are either or 1. The transformation 
will not change bits beyond these two extra places. Then hy definition 

Ei\p\iQQ)=Ei\p\ill)=Q 

and Ei acts on the remaining subspace according to the matrix formulas we have given above. This means 
that we can regard Ei diagrammatically as a cup-cap combination. E = ^ where the legs of the diagram 

correspond to the indices that can be either or 1. In this formalism ^ takes the role of and \_j the role 
of V in the decomposition E = vv^ . Each cap or cup can receive only two indices and these must be different. 



since the transformations are zero when there is a repetition of or a repetition of 1. The only further relation 
that is needed to prove from this diagrammatic point of view is that ^^l^ is an identity transformation as a 

mapping defined on a single bit. The relations EiEi+iEi = Ei follow from this. To see that this composition is 

11 



the identity, consider one of its cases as shown here: 



The fact that each cup and each cap can only 



support a zero and a one, shows that the composition will necessarily be a multiple of the identity. To see the 
details, we note that the local binary bit sequence must match the indices, which in this case is 101. This means 
that z{i + 1) = z{i) + 1 (due to the starting 1 in the local bitstring, shifting to the next node in the graph). The 

scalar contribution of ( °T J is equal to 



where the first factor under the square root comes from the cap, and the second factor under the square root 
comes from the cup. The zero in the string 101 corresponds to the z{i + 1) — 1 in the second factor. Since we 
know that z{i + 1) = z{i) + 1 it follows that S = 1. This is one of the small number of cases to check that 



proves that 



is the identity transformation. This completes the proof that $ : TL„ — > Matr{Hn,r) is a 



representation of the Temperley-Lieb algebra. 



Remark. As we have seen in the discussion above, it helps to take a diagrammatic point of view. Here we 
remark via Figure 2 on a simpler representation of the Temperley-Lieb algebra that is analogous to the AJL 
representation that we have disciisscd here. In Figure 2 we illustrate diagrams for the basic elements of the 
Temperley-Lieb algebra as we have explained them in the previous section. In order to build the elements 
Ui = Ei we make them diagrammatically - each a combination of a cup and a cap (and appropriate identity 
lines), as illustrated in this figure. In the representation illustrated in this figure, each cup and each cap is 
represented by the same matrix M and the necessary conditions for that M arc that the sum of the squares of 
its entries should be equal to d = 2cos{9) (as above), and that the matrix product of the cup and cap matrices 
should be equal to the identity matrix. This is illustrated in the middle of Figure 2 by the identity showing the 
wavy line being pulled straight to form an identity line. All the diagrams in the figure correspond to matrices, 
with the indices for matrix elements corresponding to labelled endpoints of the lines in the diagram. When two 
diagrammatic matrices are composed, an output line of one matrix is attached to an input line of the other 
matrix. Thus in these diagrams, one sums over all the possible labels of an edge that has no free ends. This 
corresponds directly to the formula for matrix multiplication where one sums the products of individual matrices 
over all double occurrences of indices. 



Remark on the Three Strand Representation. Now look at the special case of a line graph with three 
nodes and two edges: 

1 2 3. 

The only admissible binary sequences are |110) and |101), so the space corresponding to this graph is two 
dimensional, and it is acted on by Ei with z(l) = 1 in both cases (the empty walk terminates in the first node) 
and E2 with z{2) = 2 for |110) and z{2) = 2 for |101). Then we have 

E^illlO) = 0,£;i|101) = rf|101). 



E2\xyz) = \v){v\xyz) 
(xyz = 101 or 110) where v = {^/l/d, ^rf - l/d)"^. 




U'=d U 



2 M^i^M^^ = -A^-A"^ 



Figure 2. Diagrammatics for A Simple Temperley-Lieb Representation 




Figure 3. Trace Formula and Loop Count 



n 



TT 



Figure 4. Proof of the Trace Formula 



If one compares this two dimensional representation of the three strand Temperley - Lieb algebra and the 
corresponding braid group representation, with the representation Kauffman and Lomonaco use in their paper, 
it is clear that it is the same (up to the convenient replacement of A = exp(i6) hy A = iexp{i6/2)). The 
trace formula of AJL in this case is a variation of the trace formula that Kauffman and Lomonaco use. See the 
next section for a general discussion of the trace. Note that the AJL algorithm as formulated in^ does not use 
the continuous range of angles that are available to the KL algorithm, but our generalization does allow this 
continuous angular range. 



5. THE TRACE FUNCTION 

We now treat the trace function on the AJL representation from a diagrammatic point of view. Let $ : TLn — > 
Matr{Hn,r) denote the representation of the Temperley- Lieb algebra discussed in the last section. Let M = $(a) 
for any element a of the Temperley-Lieb algebra. We define a trace functional TR{M) by the formula 

Tr{M)=Y,WMk) 

k 

where Afe = sin{k9) as in the previous section, tr denotes standard matrix trace and Mk denotes the restriction 
of M to walks that end on the node k in the graph Gr- Just as in AJL, since this trace can be used to compute 
the bracket polynomial of the closure of the braid at our admissible values of A on the unit circle, it follows 
that there is a quantum algorithm for this computation by having the quantum computer evaluate the standard 
traces tr{M}~). Each Mk is a unitary matrix, and its trace can be found via the Hadamard test. 

We will prove that 

dTR{ME„) = TR{M') 

where M' is the inclusion of M in the corresponding representation of TLn+i and En denotes the matrix 
representation of the element En in TLn+i- This formula is exactly what is needed to have a Markov trace on 
the corresponding representation of the Artin braid group and hence to have a link invariant corresponding to 
the bracket model*" of the Jones polynomial. To see this from the diagrammatic point of view, examine Figure 3. 
In that figure we have shown the diagrammatic version of the formula dTR{MEn) = TR{M') at the top of the 
figure, and we have compared with the corresponding loop count formula that corresponds to bracket polynomial 
calculation. The trace we define on the representation of the Temperley-Lieb algebra will agree with the loop 
counts performed by the bracket expansion on the closure of the braid, when the TR formula is satisfied. In the 
figure LC refers to loop-count and denotes the evaluation of a collection of loops as d'^ where c is the number of 
loops in that collection. 

Finally, Figure 4 is a diagrammatic rendering of the proof of the trace formula. The expansion in the middle 
of the derivation corresponds to the evaluation of a representation of a Temperley-Lieb projector as described in 
the last section, and it takes into account the action of this projector on the bitstrings in the representation. In 
the middle of the calculation the labels and 1 remind the reader of the bitstring values for which these terms 
are non-zero. At the end of the calculation, the label k reminds the reader that these matrices are acting on 
walks that end on node k in the graph. We have left some of the details to the reader, but the main line of the 
argument is in Figure 4. The reader should note that the syntax of this use of the diagrams is explained in the 
last section. This completes the proof of the trace formula and hence the; proof that this generalization of the 
AJL algorithm to continuous angular ranges can be used as a quantum algorithm for the Jones polynonmial. 



REFERENCES 



1. D. Aharonov, V. Jones, Z. Landau, A polynomial quantum algorithm for approximating the Jones polyno- 
mial, quant-ph/0511096. 

2. D. Aharonov, I. Arad, The BQP-hardness of approximating the Jones polynomial, quant-ph/0605181 

3. R.J. Baxter. Exactly Solved Models in Statistical Mechanics. Acad. Press (1982). 

4. V.F.R. Jones, A polynomial invariant for links via von Neumann algebras. Bull. Amer. Math. Soc. 129 

(1985), 103-112. 

5. V.F.R. Jones. Hecke algebra representations of braid groups and link polynomials. Ann. of Math. 126 (1987), 
pp. 335-338. 

6. V.F.R. Jones. On knot invariants related to some statistical mechanics models. Pacific J. Math., vol. 137, 
no. 2 (1989), pp. 311-334. 

7. L.H. Kauffman, State models and the Jones polynomial. Topology 26 (1987), 395-407. 

8. L.H. Kauffman, Statistical mechanics and the Jones polynomial, AMS Contemp. Math. Series 78 (1989), 
263-297. 

9. L.H. Kauffman and S. Lins, Temperley-Lieb Recoupling Theory and Invariants of Three- Manifolds, Princeton 
University Press, Annals Studies 114 (1994). 

10. L.H. Kauffman, Knots and Physics, World Scientific Pubhshers (1991), Second Edition (1993), Third Edition 
(2002). 

11. L.H. Kauffman, Quantum computing and the Jones polynomial, math.QA/0105255, in Quantum Computa- 
tion and Information, S. Lomonaco, Jr. (ed.), AMS CONM/305, 2002, pp. 101-137. 

12. L. H. Kauffman and S. J. Lomonaco Jr. g-deformed spin networks, knot polynomials and anyonic topological 
quantum computation. J. Knot Theory Ramifications 16 (2007), no. 3, 267-332. 

13. L. H. Kauffman and S. J. Lomonaco Jr., The Fibonacci Model and the Temperley-Lieb Algebra. International 
J. Modern Phys. B, Vol. 22, No. 29 (2008), 5065-5080. 

14. L. H. Kauffman and S. J. Lomonaco, Jr., A 3-Stranded Quantum Algorithm for the Jones Polynomial, Proc. 
SPIE, vol. 6573, (2007), 65730T-1-65730T-13. http://arxiv.org/abs/0706.0020 

15. S. J. Lomonaco Jr. and L. H. Kauffman, A Quantum Manual for Computing the Jones Polynomial, Proc. 
SPIE on Quantum Information and Computation VI, Vol. 6976, (2008), pp. 69760K-1 to 69760K-4. 

16. R. Marx, A. Fahmy, L. H. Kauffman, S. J. Lomonaco Jr., A Sporl, N. Pomplun, T. SchulteHerbriiggen, J. 
M. Meyers, and S. J. Glaser, Nuclear-magnetic-resonance quantum calculations of the Jones polynomial. 
Physical Review A, (2010). 

17. R. Marx, A. Fahmy, L. H. Kauffman, S. J. Lomonaco Jr., A Sporl, N. Pomplun, T. SchulteHerbriiggen, J. 
M. Meyers, and S. J. Glaser, (In preparation.) 

18. C. Nayak, E. H. Simon, A. Stern, M. Freedman, S. Das Sarma, Non-abelian anyons and topological quantum 
computation. Rev. Modern Phys. 80 (2008), no. 3, 1083-1159. 

19. E. Witten. Quantum Field Theory and the Jones Polynomial. Comm. in Math. Phys. Vol. 121 (1989), 
351-399. 



